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Abstract 

The form factor equations are solved for an SU (N) invariant S-matrix 
under the assumption that the anti-particle is identified with the bound 
state of iV — f particles. The solution is obtained explicitly in terms of 
the nested off-shell Bethe ansatz where the contribution from each level is 
, written in terms of multiple contour integrals. 

O . PACS: ll.lO.-z; ll.lO.Kk; 11.55.Ds 

Keywords: Integrable quantum field theory, Form factors 

«D , 1 Introduction 

The Bethe ansatz , was initially formulated by Bethe 75 years ago to solve the 
^ I eigenvalue problem for the isotropic Heisenberg model. The approach has found 

^ ' applications in the context of several integrable systems in different areas, such 

as Statistical Mechanics, Quantum Field Theory, Condensed Matter Physics, 
Atom and Molecular Physics, among others. The original techniques have been 
refined into several directions: Lieb and Lininger 2^ solved the one-dimensional 
bose gas problem with (5-function potential using the Bethe ansatz. The 6- 
vertex model was solved by Lieb |31 15 with the same technique. C.N. Yang 
and CP. Yang [21 proved 'Bethe's hypothesis'^ for the ground state of the 
anisotropic Heisenberg spin chain. Due to Yang and Baxter we have 
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^Yang and Yang decided to honor Bethe's insight by calling his assumption "Bethe's hy- 
pothesis" .6 , now usually called "Bethe ansatz" . 
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the fundamental Yang-Baxter equation for the two-particle S-matrix or for the 
matrix of the Boltzmann weights, which guarantees exact integrability of the 
system. Subsequently, Faddeev and collaborators [HI E3 formulated these ideas 
in an elegant algebraic way, known as the "algebraic Bethe ansatz". Yang [7] 
and Sutherland generalized the technique of the Bethe ansatz for those 
cases where the underlying symmetry group is larger than SU{2). This method 
is now called the "nested" Bethe ansatz. This technique was applied in J2] 
to derive the spectrum of the chiral SU{N) Gross- Neveu model |13| . The 
algebraic nested Bethe ansatz was formulated in for the SU{N) and in 
jl5j for the 0{2N) symmetric case, respectively. Another generalization of the 
Bethe ansatz is the "off-shell" Bethe ansatz, which was originally formulated 
by one of the authors (H.B.) ^lElElEI to calculate correlation function in 
WZNW models (see also |2U1 \'21\ ). This version of the Bethe ansatz paves the 
way to an analysis of off-shell quantities and opens up the intriguing possibility 
to merge the Bethe ansatz and the form factor approach. In this context we 
point out that recently the form factor program has received renewed interest in 
connection with condensed matter physics |.22.. 23, 24; and atomic physics |25| . 
In particular, applications to Mott insulators and carbon nanotubes [2^1 122] 
doped two-leg ladders [27j and in the field of Bose-Einstein condensates of 
ultracold atomic and molecular gases |28[ I25j have been discussed and in some 
instances correlation functions have been computed. 

For integrable quantum field theories in 1-1-1 dimensions one of the authors 
(M.K.) et al. formulated the on-shell program |29j i.e. the exact determination of 
the scattering matrix using the Yang-Baxter equations and the off-shell program 
j3()j i.e. the exact determination of form factors which are matrix elements 
of local operators. This approach was developed further and studied in the 
context of several explicit models by Smirnov |35 who proposed the form factor 
equations (i) - (v) (see below) as extensions of similar formulae in the original 
article |[30 j . The formulae were proven by two of the authors (H.B. and M.K.) 
et al. [32] • In this article the techniques of the "off-shell" Bethe ansatz was 
used to determine the form factors for the sine-Gordon model. There, however, 
the underlying group structure is simple and there was no need to use a nested 
version of the off-shell Bethe ansatz. In the present article we will focus on the 
determination of the form factors for an SU (N) model. The procedure is similar 
as for the scaling Z{N) Ising and affine A{N — 1) Toda models |331 134j because 
the bound state structures of these models are similar. However, the algebraic 
structure of the form factors for the SU{N) model is more complicated, because 
the S-matrix possesses backward scattering. Therefore we have to apply a 
nontrivial algebraic off-shell Bethe ansatz. For > 2 we have to develop the 
nested version of this technique (see also ) . 

It is expected that the results of this paper apply to the chiral SU{N) 
Gross-Neveu model ^1 EEl EZl IHH| • In a separate article jBH] we will investigate 
these physical applications and compare our exact results with two different 
1/A^-expansions of the chiral Gross-Neveu model and [SH|- We note that 
SU{N) form factors were also calculated in [^|ini|^ using other techniques, 
see also the related paper 
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1.1 The SU{N) S-matrix 

The general solutions of the Yang-Baxter equations, unitarity and crossing re- 
lations for a U{N) invariant S-matrix have been obtained in The S-matrix 
for the scattering of two particles belonging to the vector representation of 
SU{N) can be written as 

s%{e) = SajSps h{e) + Sas^/s-r c{9) . (1) 

Unitarity reads as Si{—9)Si{9) = Ifor the S-matrix eigenvalues 

s+{e) = b{d) + c{e) , S-{e) = b{e) - c{e) . 

The amplitude Sj^{6) = a{9) is the highest weight w = (2,0, ... ,0) S-matrix 
eigenvalue for the two particle scattering. It will be essential for the Bethe 
ansatz below. 

As usual in this context we use in the notation 

^l...n g yl...n ^ yl ^ . . . ^ yn 

for a vector in a tensor product space. The vector components are denoted by 
ya _ Below we will also use co-vectors wi...n G (l/^---")^ (the dual of 

yi...n-^ with components Va- A linear operator connecting two such spaces with 

a' ...a' ~ 

matrix elements Aa\...a^ is denoted by 

where we omit the upper indices if they are obvious. All vector spaces are 
isomorphic to a space V whose basis vectors label all kinds of particles (e.g. 
V = for the vector representation of SU{N)). The vector spaces is 
associated to a rapidity variable 9i. An S-matrix such as Sij{9ij) = Sfj{9i — 9j) 
acts nontrivially only on the factors Vif^iVj and exchanges these factors. Using 
this notation, the Yang-Baxter relation writes as 

'S'l2(^12)'S'l3(^13)'S'23(^23) = 5'23(^23)5'l3(^13)5'l2(^12) (2) 

and implies here the relation between the amplitudes jlSl 

c{9) = -'^b{9) , ^ = |. 

A solution [ini Eni of all these equations writes as 

r (l - J-]T (l - ^ + 

2m) ^ \^ N 2-Ki) 

This S-matrix possesses a bound state pole in S- (9) i.e. in the anti-symmetric 
tensor channel. It is consistent with Swieca's [451 1441 138j picture that the anti- 
particle is a bound state of — 1 particles (see also jHSl El] ) • 
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For later convenience and in order to simplify the formulae we introduce 

~ _ S{9) _ 19 - P^v 

where 1 is the unit, P the permutation operator. We depict this matrix as 

5 7 



S%{9l2) = X, = Sa-,Sf3S b{9i2) + SaS^is^f c{9 




12 



and the amplitudes are explicitly 

9 — 17] 9 — 17] 

1.2 Generalized Form factors 

For a state of n particles of kind with rapidities 9i and a local operator 0{x) 
we define the form factor functions F^^ ^^{9i, . . . , 9n), or using a short hand 
notation F^{9), by 

( I 0{x) \ei,...,eny^ = e-'^'^P^+-+P-^F^{9) , for 01 > • • • > (5) 

where a = (ai, . . . , a„) and 9 = {9i, . . . , 9^)- For all other arrangements of the 
rapidities the functions F^{9) are given by analytic continuation. Note that 
the physical value of the form factor, i.e. the left hand side of ©, is given 
for ordered rapidities as indicated above and the statistics of the particles. 
The F^{9) are considered as the components of a co- vector valued function 

.n(^) e Vl...n = {V^-'^y which may be depicted as 



F^'.AO) = ■ (6) 

9i\ ■ ■ ■ 1 6*71 

Now we formulate the main properties of form factors in terms of the func- 
tions F[^ ^. They follow from general LSZ-assumptions and "maximal analyt- 
icity", which means that F^ ^{9) is a meromorphic function with respect to all 
0's and in the 'physical' strips < Im^jj < vr {9ij = 9i — 9j i < j) there are 
only poles of physical origin as for example bound state poles. The generalized 
form factor functions satisfy the following 



Form factor equations: The co- vector valued auxiliary function F^ is 
meromorphic in all variables 9i, . . . ,9n and satisfies the following relations: 

(i) The Watson's equations describe the symmetry property under the per- 
mutation of both, the variables 9i,9j and the spaces i, j = i + 1 at the 
same time 

F%...i- ■ ■ ■ ■) = ■ ■ ) ^i) 6*1, • • • ) Sij{9ij) (7) 

for all possible arrangements of the 0's. 
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(ii) The crossing relation implies a periodicity property under the cyclic per- 
mutation of the rapidity variables and spaces 



\ in, conn. 
^Pn 12. ..n 



-*'i(pi|O(0)b2,. 

= Ff'.„(^i + ivr, ^2, . . . , ^n)a?C" = F2^.„i(e2, ...Ml- ivr)C" (8) 

where takes into account the statistics of the particle a with respect 
to O. The charge conjugation matrix C^^ will be discussed below. 

(iii) There are poles determined by one-particle states in each sub-channel. In 
particular the function F^{6) has a pole at 9i2 = ivr such that 



Res Ff'.^iei,...,! 

712="r 



9n) (l — S2n • • • 5*23) . 

(9) 

(iv) If there are also bound states in the model the function F^{ff) has addi- 
tional poles. If for instance the particles 1 and 2 form a bound state (12), 
there is a pole at 6*12 = i^/, (0 < < vr) such that 

Res n{^l,92, ■ ■ ■ , On) 



^{12). ..n 
.(12) 



(^(12),...,0n)^/2r;f (10) 



where the bound state intertwiner F-^g and the values of 9i, O2, G(i2) and 
77 are given in |46[ I47j . 

(v) Naturally, since we are dealing with relativistic quantum field theories we 
finally have 

Ff^.J^l + ^, . . . , 0„ + ;,) = e^'^ i^f..„(^l, . . . , ^n) (11) 

if the local operator transforms under Lorentz transformations as O — > 
e^^O where s is the "spin" of O. 

The property (i) - (iv) may be depicted as 



(i) 

(ii) 

(iii) 
(iv) 



O 



1 cnriTi 




( ) = 


C ) 



— Res 

2l 6*12=171 



—= Res 

V2 ei2=i»7 



O 




where x denotes the statistics factor cr'-'. As was shown in the properties 
(i) - (iii) follow from general LSZ-assumptions and "maximal analyticity" . 

We will now provide a constructive and systematic way of how to solve the 
form factor equations (i) ~ (v) for the co- vector valued function F^ ^, once the 
scattering matrix is given. 
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Minimal form factor: The solutions of Watson's and the crossing equa- 
tions (i) and (ii) for two particles with no poles in the physical strip < Im 9 < n 
and at most a simple zero at 6 = are the minimal form factors. In particular 
those for highest weight states are essential for the construction of the off-shell 
Bethe ansatz. One easily finds the minimal solution of 

F{e) = a (6) F (-6) = F {2-Ki - 9) 

using (jSJ as 

oo 


It belongs to the highest weight tt; = (2, 0, . . . , 0). We define the corresponding 
'Jost-function' as for the Z{N) models by the equation 

N-2 N^l 

l[cP{e + kiri)llF{9 + kiri) = l , ri=— (13) 

fe=0 A,-0 

which is typical for models where the bound state of — 1 particles is the 
anti-particle [Hi]. The solution is 

and satisfies the relations 
</,(6») = (j){-9)a (-9) = 0((iV - l)iri - 9) 



^ -ct>i27ri - 9) = ^^^—^(t>{9 - 27Ti) . (15) 



-b{9r' ' -b{9) 

Notice that the equations (|14|) and pHj) also determine the normalization con- 
stant c in ((T^ as 

c = r-<'-'/-. (1 _ 1) exp f- f e*' f - ii^l ) m] . 

\2 2N J ^\ Jo \ tsinh^t tsinht J J 

Generalized form factors: The co-vector valued function @ for n- 
particles can be written as 

^f..n(0=i^?..n(^) n ^(%) (16) 

l<i<j<n 

where F{9) is the minimal form factor function H12|) . The K-function ^(0) 
contains the entire pole structure and its symmetry is determined by the form 
factor equations (i) and (ii) where the S-matrix is replaced by S{9) = S{9)/a{9) 

K%S ..,9i,9„...) = K%S ■■,9j,9i,...) Sij{9ij) (17) 

K?...n{Oi + iTT, 92,..., 9n)af&^ = KZ.ni{02, • • • , ^n, ^1 " ivr)C" (18) 

for all possible arrangements of the 0's. 
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1.3 Nested "off-shell" Bethe ansatz for SU{N) 



We consider a state with n particles and define as usual in the context of the 
algebraic Bethe ansatz j^l the monodromy matrix 



Tl...n,o{0_, 9q) — Sio{6io) • • • SnoiOno) 









1 


n 






(19) 



It is a matrix acting in the tensor product of the "quantum space" yi-- " = 
y^(g)- • -(gil/" and the "auxiliary space" V^. All vector spaces are isomorphic 
to a space V whose basis vectors label all kinds of particles. Here we consider 
V = as the space of the vector representation of SU{N). The Yang-Baxter 
algebra relation for the S-matrix ^ yields 



Sab{Ga 



l...ra,f)l£) ' 



' l...ra,al 



(20) 



a 






b 






b 


1 


n 




a 









b 


b 






a 




1 


n 





which implies the basic algebraic properties of the sub-matrices A, B, C, D with 
respect to the auxiliary space defined by 



Tl...nfl{§.-, z) 



2<P,P' <N. 



(21) 



We propose the following ansatz for the general form factor ^{ff) or the K- 
function defined by (|16|) in terms of a nested 'off-shell' Bethe ansatz and written 
as a multiple contour integral 



(22) 



where h{9,z) is a scalar function which depends only on the S-matrix and not 
on the specific operator 0{x) 




i=i j=i 
1 



n < 



TiZ 



Zi Zj ) 
l<j<j<m 



(23) 

(24) 



For the SU{N) S-matrix the function (j){9) is given by (|13j) with the solution 
H14|l . The integration contour Cg is depicted in Fig. ^ The constant R is 
defined hy R = <fg dz(f)(9 — z) where the integration contour is a small circle 
around z = 9 as part of Cg. 
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+ 27ri(l-i) 



92 + 2^i(l - ^) 



?i + 27ri(l-^) 



o 6*2 - 27rz;^ 



0^1 



o 6*1 - 27ri;i^ 



9n — 2711 



f 



62 - 2TTi 



6*1 - 2m 



Figure 1: The integration contour Ce- The bullets refer to poles of the integrand 
resulting from a{9i — Zj) (p{9i — zj) and the small open circles refer to poles 



originating from b{9i — Zj) and c{9i 



The dependence on the specific operator 0{x) is encoded in the scalar p- 
function p^{9,z) which is in general a simple function of e^' and e^-' (see below). 
By means of the ansatz (|16|) and (|22|) we have transformed the complicated 
form factor equations (i) - (v) (which are in general matrix equations) into 
much simpler scalar equations for the p-function (see below). The K-function 
is in general a linear combination of the fundamental building blocks |48[ \',VA\ l,S4j 
given by (|^^ - (|^H) . We consider here cases where the sum consists only of one 
term. 

If in 1)21(1 the range of /3's is non-trivial, i.e. if > 2 the Bethe ansatz 
co-vectors are of the form 



^a{0,z)=Lf3{z) $1(^,1) 



(25) 



where summation over all /3 = (/?i, . . . , /5m) with /3j > 1 is assumed. The basic 
Bethe ansatz co- vectors „ G (y^---")^ are defined as 

/3i 



(3l Pm 1 









I Oi 


9n 











2<(3i<N 
1 < a,- < iV . 



(26) 



Or 



Here the "pseudo-vacuum" is the highest weight co-vector (with weight w 
(n,0,...,0)) 

^i...n = e(l) • • • e(l) 
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where the unit vectors e(a) (a = 1, . . . ,N) correspond to the particle of type 
a which belong to the vector representation of SU{N). The pseudo- vacuum 
vector satisfies 

ni,„nBjj9,z) = 

^l...n ^l...n(^5 2) = ^l...n (^2'7) 

D^i.nM^ ^) = ft HOi - Z)^l...n ■ 

i=l 

The amplitudes of the scattering matrices are given by eqs. and ©• The 
technique of the 'nested Bethe ansatz' means that one makes for the coeffi- 
cients Lp[z) in (|^^ the analogous construction - (^1)) as for Ka{6) where 
now the indices /3 take only the values 2 < f3i < N. This nesting is repeated 
until the space of the coefficients becomes one dimensional. 

In this article we will focus on the determination of the form factors for an 
SU{N) S-matrix. The paper is organized as follows: In Section 2 we construct 
the general form factor formula for the simplest SU{2) case. In Section 3 we 
construct the general form factor formula for the SU (3) case, which is more 
complex due to the presence of the nesting procedure (more explicitly, we have 
here two levels). We extend these results in Section 4, where the general form 
factors for SU{N) are constructed and discussed in detail. 



2 SU{2) form factors 

In this section we start with the simplest case. We perform the form factor 
program for the SU{2) S-matrix. The results should apply to the well-known 
SU{2) Gross-Neveu model (TSj, investigated by Bethe ansatz methods in [i9| 
I50| . Prom the technical point of view the calculation of the form factors is very 
similar to the one done for the sine-Gordon model in [321 147j . 



S-matrix: The SU{2) S matrix^ is given by and (jSJ for N = 2. It turns out 
that the amplitudes satisfy the relations h{9) = —a{i'K — 9) and c{9) = c{i'K — 9) 
which may be written as 



= - {Sl6f, K^^ -0) + e'^e^p c(m - 9)) . 

We have introduced 



a/3 « /5 



(28) 



which are antisymmetric and £12 = e = 1 such that 



and e^f3A^€"'^= [A)^ 



tiA (29) 



^It is related to the sine-Gordon S-matrix for 13"^ = Stt by {a,b, c)^'^^^^ = (a, — fe, c)^ ''' (see 



e.g. 
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for a matrix A. Formula (P5|l may be understood as an unusual crossing relation 

(c.f. m) 

S%{0) = -C„,,5^;^(i7r - e)C^'^ (30) 




if we define the "charge conjugation matrices" as Cafs = ^afs and C"'^ = e"^. 
This means that particle 1 is the anti-particle of 2 and vice versa. 

For 9^0 and in the S-matrix yield the permutation matrix and the 
annihilation-creation operator, respectively, or in terms of 5" = S/a 

(5 7 s 'J 

a P a 13 

Ansatz for form factors: Because of (|29|) the crossing equation (|18|) writes 
in components as 

KZ...a^ (^1' ^2, . . . , en)C7^ = -K^,,,^^^, (^2, • • • , ^n, ^1 " 2m) (31) 

The form factor equation (iii) here reads as 

Res Ff „(0i, . . . , 0„) = 2i £12 Ff ^i^s, • • • , ^n) (l - c7^S2n • • • ^23) • (32) 

ai2=iTT 

Using the crossing relation ()3U() it turns out that the statistics factor in 1)18^ 
and has to satisfy {cr^) = (—1)" with the solutions 

^ -Qo foj. = nmod2 . (33) 

We define the SU{2) 'Jost-function' (j){9) and t{z) by the same equations 
as for the sine-Gordon 'Q2^ and the Z{2) models IHl] 

(l){e)F{9)F{e + i7r) = 1 (34) 
T{z)<Piz)<P{-z) = 1 (35) 

with the solution for 0(^) given by (|14|) . The functions 4>{z) and r(z) are for 
SU (2) explicitly given as 

z\„/l z\ z sinh z 



(^) = '^(--)=rfTT^)r(---^) , r{z) 



\2TTiJ \2 2m) ' ^ ' 47r3 

The scalar function h(9,z) in (|22j) encodes only data from the scattering 
matrix. The p-function p'^{9_,z) on the other hand depends on the explicit 
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nature of the local operator O. It is analytic in all variables and in order that 
the form factors satisfy (i), (ii) and (iii) of eqs. (O - © Pnm{^,z) has to satisfy 



where a is the statistics factor of ()31() and H32() . In order to simplify the notation 
we have suppressed the dependence of the p- function p*^ and the statistics factor 
cr'-' on the operator 0{x). By means of the ansatz (|16() and 1)22(1 - (|24() we have 
transformed the complicated form factor equations (i) - (iii) to the simple ones 
(i')2 - (iii')2 for the p-function. 

Theorem 1 The co-vector valued function F^{9) given by the ansatz ilh}) and 
\2}3\) satisfies the form factor equations (i), (ii) and (iii) (see ^ - ^) if the 
p-function p^{0_,z) satisfies the equations (i')2) (ii')2; (iii')2 of hHh}) and the 
normalization constants satisfies the recursion relation 



The proof of this theorem for SU{2) is similar to that for the sine-Gordon model 
in [32] and may easily obtained from the ones for SU (3) or SU (N) established 
below. 



In this section we construct the form factors for the SU{3) model, which corre- 
sponds to the simplest example where the nested Bethe ansatz technique has to 
be applied together with the off-shell Bethe ansatz. In comparison to the previ- 
ous SU (2) case, there are additional properties to be obeyed by the second-level 
Bethe ansatz function (See lemma 01 below for details). 

3.1 S-matrix 

The SU{3) S-matrix is given by (0) and (j^J for = 3. The eigenvalue 5_ 
has a pole ai 6 = irj = |i7r which means that there exist bound states of two 
fundamental particles a -|- /3 — > (pa) (with 1 < p < cr < 3) which transform as 
the anti-symmetric SU{3) tensor representation. The general bound state S- 
matrix formula |46l I47j for the scattering of a bound state with another particle 
reads in particular for the SU (3) case as 




(36) 



i'K(j){iTi)F{i'K)Nn = 2i 



'Nn-2 ■ 



3 5*^7(3) form factors 




,(pcr) 

af3 



-pip' o-') of' a' fa \cl"P'(n \ 



e 



'i2=«r/ 



(37) 
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where 6(^12) is the bound state rapidity and r/ the bound state fusion angle. The 
bound state fusion intertwiner T^^a^ is defined by 



Ressf;"'(0) = yrf'":r 



/9'a'-p(p'^) 

a/3 



13' a' 



p<Cf 



a (3 



With a convenient choice of an undetermined phase factor one obtains 



-.(per) 

Q/3 



(38) 



where F = iy^ia{iri)irj is a number. Choosing in (j3ZI) special cases for the 
external particles we calculate 

= ^(^ + - i^^) = a{iTr - 9) 

= ^(^ + l^T^HO - liTr) = b{TTi - 9) 
S^^^l^l{9) = -b{9 + \m)c{9 - \m) = c{-Ki - 9) 
which may be written as 

<m 'V) = <5j:;';'^&(- -e) + e<'-'<'U^,,c[.^ - 9) 

with the total anti-symmetric tensors Cafs-y and e""^^ (ei23 = e'^^^ = 
formulae may be again understood as a crossing relation (c.f. |45j l 



1). These 



(39) 





a (3^5 a (3^5 

or Sjjp^^ \9) = C(^ai3)S' StfS (^TT — 9)C^^'^'^'^ if we write the charge conjugation 
matrices as 



Therefore we have the relations (c.f. H29|) ') 

C„(p.)C(''-)'3 = <5f , C„(,,)^^C^(''-) = tr A 

and 



-.(po-) 



(40) 



(41) 



(42) 



These results^ are consistent with the picture that the bound state of particles 
1 and 2 is to be identified with the anti-particle of 3. For later convenience 

^The physical aspects of these facts wiU we discussed elsewhere |39| . 
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we consider the total 3-particle S-matrix in the neighborhood of its poles at 
6*12, 6*23 = irj 



Sl^^'^ (^1,^2,^3) = S^J^„{e,2)Si;. (013)5^/ (023) 

6*12 - ir] 6*23 - ir] 



e ea/37 • (,43) 



3.2 Form factors 

In order to obtain a recursion relation where only form factors for the fun- 
damental particles of type a = 1,2,3 (which transform as the SU{3) vector 
representation) are involved, we have to apply the bound state relation (iv) to 
get the anti-particle and then the creation annihilation equation (iii) 

Res F^^ .^X9i,92,03, ...,6n) = -^^2)3. n(^(i2)'^3, • • • ,0n)V2r'i2^ 

Res F^2)3...n(^(12))^3, • • • ,^'n) = 2iC(i2)3-F£.„(6'4, • • • ,6'n) 
0(12)3=^^ 

^ (1 ~ (^sSsn ■ ■ ■ S34) 

where F^g^'* is the bound state intertwiner (|nH|) (see [El El), ^(12) = f (^1 + ^2) 
is the bound state rapidity, 77 = |7r is the bound state fusion angle and C(i2)3 
is the charge conjugation matrix ()4U|). We obtain with the short notation 9 = 
(^4, • • • , ^ra) 

Res Res F^^^ ^{9) = 2iei23\^rFf J9) (l - afS^n ■ ■ ■ S34} (44) 

S23=iV 9i2=irj 

where (j42j) has been used. 



Ansatz for form factors: We propose that the n-particle form factor of 
an operator 0{x) is given by the same formula (|16|) as for SU{2) where the 
form factor equations (i) and (ii) for K- function write again as H17|) and (|18() . 
Consistency of (|44|) , (|18|) and the crossing relation (|39jl means that the statistics 
factors are of the form cj^ = a'-'{ra) if the particle of type a belongs to a SU{S) 
representation of rank = 1,2 and 

a^{r) = e^'^i^^o for Qo = nmod3 (45) 

as an extension of (|33|) . As for SU{2), we propose again for the K-function 
the ansatz in form of the integral representation ()22() with H23|). However, the 
Bethe ansatz co-vector is here for SU{3) of the form 

^a{^,z)=Lp{z)^^{9,z) (46) 

where the basic Bethe states ^^{9, z) are given by (PH)) and the indices run 
over 1 < < 3 and the Pi over 2 < /3j < 3. For the function Lp{z) one makes 
an analogous ansatz (|22() as for Ka{ff) where the indices run here over a set 
with one element less. By this procedure one obtains the nested Bethe ansatz. 
The next level function Lp{z) is assumed to satisfy 
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(i) (i) 

• • 1 1 • • • ) — L ji Jy. . . J Zj , Zij . . . ) Sij{^Zij) (4"^) 

(ii) (i) 

C^^Li..,m(zi + ivr, Z2, . . . , Zm) = L2...ml{z2, ■ ■ ■ , Zm, Zi - Z7r)C^^ (48) 

(iii) ^^) there is a pole at 2:12 = ii] such that 

m 

Res = ci 170(2,2) Res Sf.{zi2)Lg{z) (49) 



2 = 3 



with ci = (f){iri). 

These properties of the next level Bethe ansatz function Lp{z) are discussed in 
lemma |31 

The minimal two particle form factor function 

00 



belongs to the highest weight w = (2,0,0). The SU{3) and the Z{3) model 
[SSI ISl] possess the same bound state structure, namely the anti-particle is to 
be identified with the bound state of two particles. Therefore we define the 
SU{3) 'Jost-function' (j){z) by the same equation as for the Z{3) model 

(l)iO)(l){0 + ir])F{e)F{0 + i7])F{e + 2ir]) = 1 (51) 

with the solution 

' \2m) \3 2m 

These equations determine also the constant c in (|5()|) . The function r(z) is 
again defined by PH) . 

The function h{9,z) is scalar and encodes only data from the scattering 
matrix. The p-function p'^{0,z) on the other hand depends on the explicit 
nature of the local operator O. It is analytic in all variables and, in order that 
the form factors satisfy (i) - (iii) of eqs. © - @, we assume 



(i')3 p{§.,z) is symmetric under 9i ^ 6 

.../^ / ap(^i + 27ri,02,...,l) = (-l)X^i,e2,...,l) 

1 p(^,zi + 27ri,Z2,...) = (-l)>(^,^i,^2,...) (52) 



(iii')3 for 6*12 = 6*23 = irj 



p{0;e2,e^,z) = {-irp{e,z) 

p{O;Oi,92,z)=ap{0,z) 



with 9 = (6*4, . . . , 9n), z = (23, . . . , Zm)- In (ii')3 ^^'^ (iiiOs ^ is the statistics 
factor of the operator O with respect to the fundamental particle belonging to 
the vector representation of SU{3). 
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Theorem 2 Let the co-vector valued function F^{6_) be defined by the ansatz 
ilb]} . \2^) . \2'J\) and \4(^ - Let the p-function satisfy (i')3, (ii')3 ^'^'^ (iiiOs ^/ 
iTTl)) and let the function Li3{z) satisfy (ii)^^^ and (iii)^^) of - ^^ >. 

Let the normalization constants satisfy the recursion relation 



2{ir]f4>'^{iri)4>{2ir])F'^{ir])F{2ii])Nn = 2i^/2r7V„_3 . 



(53) 



Then the function F^{0) satisfies the form factor equations (i), (ii) and (iii) 
(see - In particular i44\ ) is satisfied. 

Proof. Property (i) in the form of H17() fohows directly from (i')3, the 
Yang-Baxter equations and the action of the S-matrix on the pseudo-ground 
state Vt 

^...ji... C^_ji^^^[. . . 6j , 6i . . .) SijiOij) = SijiOij) C_^ij^^\. . . 9i, 6j . . 

— ^...ij... C ...ij ..X' ' '^i^^j' ' •) ■ 

because 311(9) = Sll{9)/aie) = 1 and F{9) = F{-9)a{9). 

Using (i) and (|4H) the property (ii) in the form of (|18j) may be rewritten as 
a matrix difference equation EH ESI 



Ki,„n{9:)a = Ki,„ni9)Qi...n{ 



(54) 



where 9' = {9i + 2Tri, 92 ■ ■ ■ , 9n) and cr is a statistics factor. The matrix Q{&) 
1) is a special case (for i = 1) of the trace 



Qi...n{9,i) = tro Tq 



,l...n,0 



OL,; 



OLi 



Ctr. 



of a modified monodromy matrix 



!Q,l...n,0(£, «j 



S 



■ PiO ■ ■ ■ Sno{9n 



where P is the permutation matrix. For the special case i = 1 the matrix 
Qi^i 1) = Q{^) may be written as a trace of the ordinary monodromy matrix 
H19|) over the auxiliary space for the specific value of the spectral parameter 
^0 = ^1 

Qi...n(^) =trori...„,o(^,^i) (55) 
because 5(0) = P. The difference equation may be depicted as 



K 



K 
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where we use the rule that the rapidity of a hne changes by 27ri if the hne bends 
by 360*^ in the positive sense. 

In the following we will suppress the indices 1 . . . n. The Yang-Baxter rela- 
tions H2U() imply the typical commutation rules for the matrices A,C,D defined 
in eq. ^ 



-A{9, e)c^{e, 



b{9-z) b{e-z) 

^—-sj^^ iz - 9)w;{e, 9)cP\e, z) 

o{z — Uj ' 
b{z - 9) 



A{9,z)C^{9,9) (56) 



D^i9,z)C'^ (9,9) 



where /3, 7, 7', 7" S {2,3}. In addition there are the Zapletal commutation 
rules ESI ES] where also the matrices Aq , Cq , Dq defined by 



are involved |51j 

&i9,z)AQ{9) 
&{9_,z)DqX{9) 



1 



Aq{9_) B^q{9) 



-Aq{9)&{9!,z 



^^'^ iz 



m - < 
1 

~i^^m9,z)ct{9). 
h{z-9i) ' 



b{9[ - z 
[)DqX{9)&'{9_!,z 



mz)C^oi9_) (57) 



(58) 



Note that we assign to the auxiliary space of Tq{0) corresponding to the hori- 
zontal line the spectral parameter 9i on the right hand side and 9'^ = 9i + 27rz 
on the left hand side. 

We are now going to prove (|54|) in the form 



K{&)\Aq{0) + Y,D^ 

/3=2 



K{9) Q{9) = K{9!)a 



(59) 



where K{9) is a co- vector valued function as given by eq. (|22|) and the Bethe 
ansatz state H25() . To analyze the left hand side of eq. H59() we proceed as follows: 
We apply the trace of Tq which is Aq + X]/3=2 ^Qp co-vector K{9) 

/3l /3m 1 1 1 



■■■&^{9_,z,)fQX 















Zl 










02 


On 



1 

■7 
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In the contribution from Aq{6) which means 7' = 7 = 1 one may use Yang- 
Baxter relations to observe that only the amplitudes Sll{9i — Zj) = 1 appear in 
the S-matrices S{6i — Zj) which are constituents of the C-operators. Therefore 
we may shift all Zj-integration contours Cq to Cq' without changing the values 
of the integrals, because there are no singularities inside CqU —Cqi (c.f. Fig. P). 
Using a short notation we have 



K{e)AQ{9)= / dzh{e,z)p{e,z)^{e,z)AQ{e) 

JCgl 

(with J^^ dz_= J^^ ^ ' ' ' Ice ^^)- proceed as usual in the algebraic 

Bethe ansatz and push the Aq{6) and Dq{9) through all the C-operators using 
the commutation rules H57|) and H58|) and obtain 

m ^ 

&-{9, z^)... &^{9, zi) Aq{9) = n r-, r Mi) 

j=i ^j) 

X &^-{9:, z^)--- cf'{9:, zi) +Y,uwA, 



III 

&-{9, z^)... C^HO, z,) Dq^^{9) = n —T^'^ 0[)Dq^;{9) 

X C"^'" {9_,Zm)--- {9, Zi)+J2UWD. 

The "wanted terms" written out explicitly originate from the first term in the 
commutations rules (|57|1 : all other contributions yield the so-called "unwanted 
terms" . The next level monodromy matrix is 



f^'^ Tp U, 0) = (5io(zi -9)--- ~Sm0{Zm - 9)) 



where the /3's and also the internal summation indices take the values 2, 3. If 
we insert these equations into the representation (|22|) of K{ff) we first find that 
the wanted contribution from Aq already gives the result we are looking for. 
Secondly the wanted contribution from Dq applied to Q gives zero. Thirdly 
the unwanted contributions from Aq and Dq cancel after integration over the 
Zj. All these three facts can be seen as follows. We have 



^Aq{9) = ^, VtDQl{9) = 5''^^b{9i-9i)n = {) (60) 



i=l 

which follow from eq. H27|) . Therefore the wanted term from Aq yields 



iw r ___ 1 _ 
m)AQ{9) = dz'[\^^ -h{9,z)p{9,z)^{9^,z) 



a / dzh{&,z)p{9^,z)-^{9^,z) 

JCgl 
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It has been used that the 'shift relation' of the function (piO) in H15|l and (ii') of 
H52|) for the function p imply that 



n 



hii,z)p{e,z) = (jh{e;,z)p{e;,z) 



The wanted contribution from Dq vanish, since 6(0) = 0. Therefore it remains 
to prove that the unwanted terms cancel. The commutation relations H57() 
and (|58)) imply that the unwanted terms are proportional to a product of C- 
operators, where one C{9,Zj) is replaced by Cq{9). Because of the symmetry 
(i)^ of Lp{z) it is sufficient to consider only the unwanted terms for j = 1 
which are denoted by uw\ and uwj^. They originate from the second term in 
the commutation rules (|57|) when Aq (6) is commuted with C {9, z\ ) . Then the 
resulting A{9_^z\) pushed through the other C's. Taking into account only the 
first terms in (|56j) we arrive at 



uw\{z] 



aA{9,zi)C^-^{9,Zm)...C^C} 



Using (|27|1 and Yang-Baxter relations (always taking into account the ^SU{N) 
ice rule' which means 'color conservation') and d{9) = 1 one obtains 



nA{9,zi)C^"^i9,Zm)...Cg'{9) 



1 1 



Zl 



Z2 



1 1 



Ql 02 



dr. 




Z2 



a\ a2 



Or, 



where ^^{9_,z) does not depend on ai, l3i,9i, zi. Similarly, we obtain the first 



unwanted contribution from Dq^ as 



C(Zl 



b{zi 



1 



'i) J=2 ^(^i 



Zl, 



{fW'^%U,z,)) 

X nD^^"{9, zi)Cf^'^{9,z„ 



c5m 
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which may be depicted as 



Zm,) ■ ■ ■ Cr)^ { 



P2 Pm ■ 


[ 1 


1 
















Zm 






















02 


On 



1 
1 
1 



a\a2 



Again (|27|) and Yang-Baxter relations yield 

^ of {e, zi)&'^ {e, z^)--- c^q' (0)] ^ = 5f J] m - Zl) 4>| {e_, 1) 

- i=l 

while assumption (ii)'^^^ for Lp{z) means 

Lp_U)Q%{z) = Lp_{^) 

where analogously to eq. ([55]) Q'^^'^^,{z) = T^^'^^^r^{z, zi) is the next level Q- 
matrix and z' = (z'l = zi + 2Tri, . . . , Zm)- Therefore we finally obtain 

"I uw 

L J a 

= - [ dzh{e,z)p{e,z) f^}~''\ f[ ^ L^uX^f (1,1) 

JCg, b{e[- Zl) fj-^b{zi- Zj) - 



uw 
a 



dzh{9_,d)p{0.,z)i^^ 



n 



h{zi - Oi) fj^ b{z[ -Zj) - 



It has been used that the 'shift relation' of the function ^{0) in (|15)1 and (ii') of 
H52|) for the function p imply 

n^(^^ - n I— — T H0,i)p{e,z) = n j-r—: m^)pio,z:) . 

i=l j=2 '^y^j j=2 "v^i ^j) 

For the Dg-term we rewrite the zi-integral by replacing zi ^ zi — 2iTi (such 
that z'l — > Zl and Cg ^ Cq + 27ri) and obtain for the sum of the unwanted term 
from Aq and Dq (using c{z)/b{z) = —c{—z)/b{—z)) 



Km [AQie) + DQie) 

c{zi - 



dz 



Ce+2iTi JC, 



dzi / dz 



h{zi 



^mi}p''{e,z)Lp_{z)\{^ 



i=2 



b{zi - Zj) 
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Taking into account the pole structure of ^^^^—^h{9,z} we may replace for the 
An -term 



dz ■ 



-2TTi Res + / dzi / dz_- 

^l=^i JCg I JCg, 



where again ( /c , ~ j • • • = and for the Dg-term 



/ dzi I dz_ 

JCg+2ni JCg 



-27rz Res + / dzi / dz--- 

^l=e[ Jcg / Jcg 



which proves the cancellation of the unwanted terms. 

The proof of (iii) is similar to that for the Z{3) model in j^. We use the 
short notations 



6 = (01, ... , 9n), = {6i,92,0'i), 6 = (6*4, . . . , 0„), 
a = (ai, . . . , On), a = (ai, 02, as), « = (04, • • • , 

Z={zi,..., Zm), Z = {zi,Z2), Z= {Z3,..., Zm), 

and prove (jlH) which may be depicted as 



Res Res 

923=iri di2=iv 



o 



2i\f2T 



O 



O 



n r 



We will show that the K-function given by the integral representation (|22|1 with 
1^5]) satisfies 



Res Res 



n 3 
i=4j=2 

Co = 2i\/2rF~2(ir/)F"^(2i?7) . 



£123 i^4...n® (1 - aa^sn . . . 534) (61) 



This is equivalent to (|44|) because of the ansatz for Fi n{6) and the relation 
of F{z) and (/)(z) given by (|51() . The residues of -firi...n(^) consists of three terms 



Res Res Ki.„„(0 = iJ^., + R^',^ + • 

6i2=iV 62:i=iV 

This is because each pair of the z integration contours will be "pinched" at 
three points. Due to symmetry it is sufficient to determine the contribution 
from the zi- and the Z2-iiitegrations and multiply the result by m{m — 1). The 
pinching points are 

(1) Zi = 02, Z2 = 03, 

(2) Zi = 01, Z2 = 02, 



?{2) 



>(3) 
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(3) zi 



The contribution of (1) is given by zi-and 2;2-integrations ifg^ dzi <fg^ dz2 ■ ■ ■ 
along small circles around zi = 62 and Z2 = 03 (see figure P). The S-matrices 
3(02 — zi) and 5(^3 — Z2) yield the permutation operator 5(0) = P 



P1P2I33 /5m 1 1 1 1 



n c"- (0, zm)--- c^^{e, 63)0^' (0, 02, 



Oir. 



Ql 0203 04 

m 

= ~s^JX^0i3ysi\il{0i2) n b{0i - zj) [n&-{0, Zm)--- c^Kl zz. 

We have used the fact that because of the SU (N) ice rule only the amplitude 
b contributes to the S-matrices ^(^i — Zj) and only a = 1 to the S-matrices 
S{02 — Zj),S{03 — Zj),S{0i — zi),S{0i — Z2) after having applied Yang-Baxter 
relations. Further we use that for Z12 023 iij by assumption (iii)^^) (c.f. 
(jH) 



L^{z) ^ Cl'[[4>{Zi2)Sf^p^{zi2)L^{z) 
i=3 

(f>{ir]) 



ci 



^12 



23) 



and that because of 



Res Res </)(0i2)<A(^i3)5;^f,,,3(^i, ^2, ^3) = </'(ir?)</'(2ir/)2 (ir/)^ e,,„,,3 . 

We combine this with the scalar functions h and p and after having performed 
the remaining Zj-integrations we obtain 

n 3 

4^)= conn 

i=4j=2 



because of b{0i-z)4>{0i-z)4>i02-z)4>{z-03)4){03-z)T{02-z)T{03-z) = -1, the 
relation (iii')3 of H52|) for the p-functions p{0_, 02, ^3,^)|6'i2=623=«'7 ~ 

i-l)"'p{0,z) 

and the relation (|53|) for the normalization constants NnCi2{iri)'^ cl){iri)(j){2ir]) = 

Nn-3C0. 

The remaining contribution to (|61l) is due to R\ + Rl It is convenient 
to shift the particle with momentum ^3 to the right by applying S-matrices 
using (i). Then we have to prove that 



n 3 



Rl.L + Ri.L) ^43... Sn3 + CO n n "^(^^i) ^123 i^4...n(0CT3 

i=4i=2 



0. 



(62) 
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Note that because of (i) (c.f. ((T7|) ) 

n 

Kl...n{e)SiZ . . . 5n3 = n 

a{9i-i)Ki„A...n'i{0i, . . . ,04, . . . ,0n,O3) . 

4=4 

Due to Lemma El in appendix^ it is sufficient to prove equation (|62|) only 
for 03 = 1 since the left hand side of this equation is a highest weight co- vector. 
This is because Bethe ansatz vectors, in particular also off-shell Bethe ansatz 
vectors are of highest weight (see 0,;). Therefore we consider this equation 
for the components with 03 = 1. The contribution of pinching (2) is given by 
the zi-, 2;2-integrations along the small circles around zi = 61 and Z2 = 62 (see 
again figure P). Now S{6i — zi) and S{92 — Z2) yield permutation operators P 
and the co-vector part of this contribution for 03 = 1 is 

{nc''"^{e,zm) ■ ■ ■ &Hd,02)c^H0,0i)P3{i))^ 

/3l/?2/9m 111 1 




ai 0.2 a4 a3 = 1 



= StXl {^C^-'iO, Zrr.)--- C^'il ^3))/a3 (63) 

where -P3(l) projects onto the components with 03 = 1. We have used the 
fact that because of the SU{N) ice rule the amplitude a only contributes to 
the S-matrices S{9i — Zj), 5(6*2 — Zj), 5(^3 — Zj), S{9i — zi) after having applied 
Yang-Baxter relations. One derives the formula 

Res Res 0(03i)<^(^32)532 (^12) = -<^(«r?)0(2i??)2(ir?)2e^^^^i 

d2'i=iri 9i2=vri 

using e''' = 1, Sfl\^^{Bs + 2^,61,62) = 6(^31 + 27ri)6(032 + 27ri)532^^(0i2) 
(because of the SU{3) ice rule) and 0((9) = -b{9 + 2TTi)4>{e + 27ri) (see We 
combine this and H63|) with the scalar functions h and p taking into account the 
property (iii)^^) of (|4,'-{)1 and after having performed the remaining Zj-integrations 
we obtain 

n 3 

i?S'^.„543 . . . 5„3^'3(1) = -CO n n "^(^'j) ^123 i^4...n(0a3P3(l) • (64) 

i=4i=2 

We have used the following equations a{6i3)(l){9ii) = (j){9i3) (because of (|T5|) ). 
the definition (pi)) of t{z) which implies (/>(z — 92)(j){9i — z)4){92 — z)(p{93 — 
z)t{9i — z)t{92 — z) = 1, the second relation (iii')3 of for the p-functions 
p{9_,9i,92, z)\ei2=d23=iv ~ '^P(^)l) the relation for the normalization 
constants NnCi(t>{ir])(f>{2ir])2{ir]y = Nn~3Co. 
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The contribution of pinching (3) is given by the zi-, 2:2-iiitegrations along 
the small circles around zi = 92 — ir] and Z2 = 9-^ — irj, (see again figure Now 



S{93 - Z2) yields S\ 
state 



(57 , 



■Z2) 



-.(57 -pCpcr) 



and the residue of the Bethe ansatz 



/3l/?2/5r> 



Res {^C{9_,Zm)---C{9,z^)Pr,(\) 



1 1 1 



1 

a\ a2 as an = 



vanishes for 03 = 1 because T^^'' is antisymmetric with respect to a, (3. There- 
fore equation (|62() is proved for 03 = 1 and because of Lemma El also in general. 



The higher level Bethe ansatz 

Lemma 3 Let the constant R and the contour Cz he defined as in the context 
of h2^). Then the higher level function 

= hl ^ • • 7 ^ ^) ^) 

J Cz Cz 

^f{z,u) = (7(1) (1, ufe) • • • C(i)(z,ni))^ 

m k 

h{z,U} =^'Y\4>{Zi - Uj) JJ T{Ui-Uj). 
i=l j=l ^<i<j<k 

satisfies (i)(i) - (iii)(i) of the equations - if 
{i')^^^ '— ) symmetric under zi ^ Zj 

{u')f^ f + 2vri,Z2,...,u) = (-l)Vli^(zi,Z2,...,w) (65) 

Proof. The proofs of the second level equations (i)^^) and (ii)*-^) are similar 
to the ones for the first level. For the proof of (iii)''^) we observe that for 
Z12 ir] there is a pinching at Uj = Z2 which means that in ^^0\z_,u) we may 
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replace S{z2 — ui) P and we have to consider (for 2:12 itj and ui = Z2) 

2 2 2 2 



/3 3 
3 











Uk 










U2 






m = Z2 




Zl 




Z3 





Pi P2 Ps 

k 

i=2 



with P = {P3, . . . , Pm), z = (2:3, . . . , Zm), u = {u2, • • • , Uk)- Therefore because 
of b{zi — u)(f){zi — u)(p{z2 — u)t{z2 — u) = —1 and (iii')3"'^'* for Z12 ir] 



Lp{z) p. ^{z^2)\{kz^2ySf^pM^)Lfz) 
i=3 

which proves the claim (iii)^-'^) of (|49() . ■ 

4 SU{N) form factors 

In this section we perform the form factor program for the general SU{N) S- 
matrix. For this purpose, we extend the procedures of the previous section, 
i.e., the nested Bethe ansatz method now with — 1 levels combined with the 
off-shell Bethe ansatz. Applications of the results to the SU{N) Gross-Neveu 
model 55 "^ill be investigated in a separate article [HH] . 



4.1 S-matrix 

The SU{N) S-matrix is given by and ©. Again, the eigenvalue S-{6) 
has a pole at 6 = irj = j^iir which means that there exist bound states of r 
fundamental particles ai + ■ ■ ■ + Ur — > {pi . . . pr) (with pi < ■ ■ ■ < pr) which 
transform as the anti-symmetric SU{N) tensor representation of rank r, (0 < 
r < N). The masses of the bound states satisfy rrir = rriN-r which suggests 
Swieca's |45| I44| I38j picture that the antiparticle of a particle of rank r is to be 
identified with the particle of rank N — r (see also !33, 34,). 

Iterating — 1 times the general bound state formula one obtains for the 
scattering of a bound state {P1P2 ■ ■ ■ Pn-i) with another particle 5 (analogously 
to dSZl) for the SU{?>) case) 



f(5'(7i72---7iV-i)//i\p(/3i/32.../3jv-i) 
{0i02-I3n-i)& ^ ' «i"2...a]v-i 
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with the total bound state fusion intertwiner 



p{/3i/32.../3iv-i) _ ■p(/3i/32.../3]v-i) t.(/3i/32/33) t.(/3i/32) 

i aia2...«jv-i - i (/3i/32.../3iv-2)"iV-i ' ' (/3i/32)«3 "i"2 • 

Taking special cases for the external particles we obtain 

Sii[i^.;N-[)N(^) = ^(^ + i7r-iv)... h{9 -m + ir,) = - 9) 



S. 



tt'N-i)N-m = + ivr - zry) . . . a{e - ivr + ir/) = {-lf-^h{7,t - 9) 
^S.'.!'ivH~^V) = -KO + iTT-ir,)... c{e - ivr + ir,) = {-l)^-^c{Tri - 9) . 
which may be written as 

^/3jv(/3i---/3jv_i) ^^-j 



(ai...ojv-i)ojv 



"ai...ajv_iajv 



civrz 



with the total anti-symmetric tensors eQi...oj^ and g"i---"iv (g;^ jy = = \y 
These results may be interpreted as an unusual crossing relation 

e/3i/32.../3iv-i/3 'S'aj^j {6l) ■ ■ ■ S^^~^^^''~'' {6n-i) = (-I)^""^eaia2...«iv-17 ^Ip (^^ " ^) 

(66) 




(-1) 



N-l 




7 



5' 



with 0j = 9 + in — jirj if we write the charge conjugation matrices as 

C{ai...aisi-i)aisi Cai(a2---ajv) €ai...ajv ) 
Q(ai...a]V-.i)a]v _ ^ai(a2---«jv) _ gOi-.-ajv 

Therefore we have the relations (c.f. (|29|) and (|41|) ) 

/-I /-i(ai...ojv_i)/3 _ e/3 f-1 /.ap/3(ai...ojv_i) _ ^.Af-l 

'-^a(ai...a]V_i)'-" — ) '-^a(ai...a]V_i)^/3'-" — \ ^) ^ 



f~i -p(/3i---/3]v-i) 
ana »^(/3^.../3^_^)^i Q,i...a^_-^ 



':ai...a]v-l7-' 



(67) 



where the constant T is 



V 27r V iV 



These results are consistent with the picture (c.f. [3H1) that the bound state 
of — 1 particles of rank 1 is to be identified with the anti-particle of rank 
1. As already remarked, the physical aspects of these facts will be discussed 
elsewhere jHS]- 
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For later convenience we consider as a generalization of H43() the total A^- 
particle S-matrix (consisting of N{N — l)/2 factors) in terms of S 



S12...N = ySuSis . . . SiNj (^5*23 . . . S2Nj ■ ■ ■ Sn~in (68) 
in the limit Ojj+i ii] {j = 1, . . . N — 1) which behaves as 

t>i2 - IT] t/N-iN - iri 

The algebraic structure of this relation follows because one can use the Yang- 
Baxter equations to shift in (|68|) any factor Sa+i ~ (1 — P)jj_|_i to the right or 
the left. Therefore the expression is totally anti-symmetric with respect to the 
and the f3i. The factor follows from Q. 

4.2 The general form factor formula 

In order to obtain a recursion relation where only form factors for the fundamen- 
tal particles of type a (which transform as the SU {N) vector representation) 
are involved, we have to apply iteratively the bound state relation (iv) to get 
the {N — l)-bound state which is to be identified with the anti-particle 

Res ... Res F^^{e) 



9i2=iri eN-2N-i=iri 



P(l...N-l)N...ni^{l-N-l)idN, • • • , 0n)V2 r[^";^_^^^ 



and finally the annihilation residue equation (iii) 

Res Fn...N-l)N...ni^il...N-l),GN,- ■ ■ ,0n) 

= 2«C(i jv_i)Ar-F^+i, „(0Ar+l, . . . ,9n) - O-'^Snu ■ ■ ■ SnN+i) ■ 

Similar as for = 3 we obtain with (|E7I) 
Res . ... Res Fi^3...n(^i, • • • , ^n) 

9n-in=iv 6i2=ir] 

= 2iei...ArA/2 TF^_^_i_^{9n^i, . . . ,6n) {1 — cr'^SNn ■ ■ ■ Snn+i) ■ (70) 

Ansatz for form factors: We propose the n-particle form factors of an 
operator 0{x) as given by the same formula H16() as for SU{2) and SU{3) in 
terms of the K-function and the minimal form factor function F (9) given by 
(|12|) which belongs to the highest weight w = (2,0, ... ,0). The form factor 
equations (i) and (ii) for K-function write again as (|17j) and (|18|) . Consistency 
of H70|) . H18|) and the crossing relation (|66() means that the statistics factors 
are of the form cj^ = a'~^{ra) if the particle of type a belongs to an SU{N) 
representation of rank = 1, . . . , N — 1 and 

a'^{r) = e*^(i-i/iV>«o for Qo = nmodiV (71) 
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as an extension of and For the K-function K^{9) we propose again 

the ansatz in form of the integral representation 1)221) with (|23() and (|24|) . The 
Bethe ansatz co- vector ^z) is again of the form ()46)) and for the function 

Lf^[z) one makes again an analogous ansatz as for the K-function 1)22(1 where 
the indices run over sets with one element less. For the SU (N) case we have to 
iterate this N — 2 times 

4' = ^ / ^ • • • / ^ M^, U, I.) U, I.) (72) 

= n n ^^^^ ~ ^j') n '^^^^^ ~ 

i=l j=l ^<i<j<k 

for / = 1, . . . , iV - 2, ai = /,..., iV, /5i = / + 1, . . . , iV, = By this 

procedure one obtains the nested Bethe ansatz^. 

As for the Z{N) and A{N - 1) models [HSl El the ' Jost-function' (j) {9) = 
<p{9)/a (9) = (f){—9) is a solution of the equation 

N-2 N-1 

'[l(t>i9 + kir]) '[[F{9 + kirj) = 1 (73) 

fc=0 fc=0 

which is typical for models where the bound state of — 1 particles is the 
anti-particle. The solution is 

' \2-Ki) V N 2iTi J 
and again we define r(z) = l/{(p{z)(l){—z)). 

The higher level Bethe ansatz: In order that the form factor F*-^ (9) given 
by (|16() and ()72() satisfies the form factor equations (i) - (iii) the higher level 
L-functions L^'^^ ^^{zi, Z2, ■ ■ ■ , Zm), {I < Pi ^ N, m = ni) have to satisfy: 

(i)(0 Watson's equations 

L {. ■ ■ , Zi, Zj, . . . ) = L [. . . , Zj, Zi, . . . )Sij[zij) , 

the crossing relation 

-^i3i,/32,...,/3m(^l' ^2, ...,Zm)= -C'^2,...,/3™,/3i ^^^^ , ■ ■ ■ , Zm, Zi - 2m) 

and 

''in |85) the nested ofF-shell Bethe ansatz was formulated in terms of "Jackson-type inte- 
grals" instead of contour integrals. 
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the function lI^^ {z) has to possess simple poles at 2:12, • • • , zn^i^i^n-i = 
irj and it has to factorize in the neighborhood of these poles as 

m N~l 

Lf{z)^ci n X{kz^,)~Sf■'^\i)L|{z) (74) 

i=N-l+l j=2 - 



L(0 



. l+l 

L(1) j mni N-l 

~TT ^ n n '^(^^i) \i) 

ft i=N~l+l j=2 I • • • I 

where we have used the short notations /? = . . . , Pn-i)i 

§_ = {(3n-i+i, ■ ■ ■,(3m) and z = {zi,. . .,zn-i), z = {z^^i+i, . . .,Zm)- The 

S-matrix S^"'^~^^{z) describes the total scattering of the — I particles 

with rapidities z for the initial quantum numbers /? and the final config- 
uration of quantum numbers (A^, ...,/ + 1). 

The constants satisfy the recursion relation 

N-l-l 

Q = Q+i n '^(•?*^) C^^) 

with the solution 

ci = 4>^-'^{iri)4>^-^{2iri) ■ ■ ■ 4>{{N - 2)ir]) (76) 
if cn-i = 1 (see appendix lU)) . 



Lemma 4 The higher level function Lp\z_) satisfies (i)^'-* - (iii)^'^, if the high 
level p-function in the integral representation \7}3{j satisfies 



(i')(0 


p«(. 


(ii')(0 


r 1 
) 

I p[ 


(iii')(') 




with m = 


M - 



mfe(^i +2vr«,22,...,2i) = (-l)V™fc(^l'^2,---,n) 
'^2k{z,ui + 2m,U2,...) = (-l)™pJJ^\(i,'Ui,'U2,---) 

:2,....^_.n)|.,,=...=..__^_=., = (-l)'=-^-'+^P« (i,^) 

(77) 



This lemma is proved in appendix [0 



The p-function: In order that the form factor F^ (9) given by 1)16^ and 1)22^ 
- (|24j) satisfies the form factor equations (i) - (iii) the p-function p^{9,z) in 
which depends on the explicit nature of the local operator O is assumed 
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to satisfy 

(i') PnmiS,z) is symmetric under 9i ^ 9j 

apnm{01 + 27Ti,62, . . . , z) = (- 1)'"+^- , ^2 , ■ ■ ■ , z) 
Pnm{0, Zl + 2m, Z2,...) = (-l)>nm(^, Zi,Z2,...) 
(iii') Pnm(^, 6*2, • • • , ^7V,l)|0i2=---=6»]v_ijv=ir? = ( — 1)™""^^"^ Pn-Afm-Af+1 l) 
PnmiS., ^Af-lil)|ei2=---=6»]v-ijv=j»7 = ^ Pn-Afm-Ar+l 1) 

(78) 

where ^= (6*1, . . . , 6'„), ^= (6'iv+i, . . . , 6'n), 1= and 
z = [zjy, . . . , Zm)- In order to simplify the notation, we have in these equations 
suppressed the dependence of the p-function and the statistics factor on 
the operator 0{x). 

Theorem 5 The co-vector valued function Fa{ff) given by the ansatz Ub]) and 
the integral representation \2S^) satisfies the form factor equations (i), (ii) and 
(iii) of ^ - 13), in particular \7(J{ ) if 

1. Lp{z) satisfies the equation {})^^\ (ii)''"^^; (iii)^^) and \7b\ l of lemma^ 

2. p^{0,z_) satisfies the equation (i'), (ii') and (iii') of ( |7($D and 

3. the normalization constants satisfy 

II [4>{3iv)F{3iv)) Nn = 2iV2 ~ TN^-n ■ 

(79) 

The proof of this theorem can be found in appendix IbI 
4.3 Examples 

To illustrate our general results we present some simple examples. 

The energy momentum tensor: For the local operator 0{x) = Tf"^{x) 
(where p,a = i: denote the light cone components) the p-function is, as for the 
sine-Gordon model in i47i 



p^''^ie,z) = Y,e^'^^e'^^^ 



i=l i=l 



For the n = N particle form factor there are ni = N — I integrations in the Z-th 
level of the off-shell Bethe ansatz. The SU{N) weights are (see 1351 ) 

w = {n- ni,ni - n2,. . . , nN-2 - nN-i,nN-i) = (1, 1, . . . , 1, 1) . 

We calculate the form factor of the particle a and the bound state (/?) = 
(/?!, . . . , Pn-i) of N — 1 particles. In each level all integrations up to one may 
be performed iteratively using the bound state relation (iv) (similar as in the 
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proof of theorem |SJ . Then all remaining integrations in the higher levels can 
be done by means of the formula 

/ dsr{a+s)T{b+s)T{c-s)T{d-s) = 2^^^ ') ' ^ . 

J-ioo T{a + b + c + d) 

The result for the form factor of the particle a and the bound state (/?) writes 
as 

F^[;^{e^,62) = Kl'^;){9^,e2)G{6^2) 

Kll^iOM = Nr {e^'^+e^'^) j^^ ^{0, - z)e'^^m " (80) 

where the summation is over 7 and 6 > 1 and G{9) is the minimal form factor 
function of two particles of rank r = 1 and r = N — 1. The functions G{9) and 
L{9) are given by 

G{i7r - 9)F{9)(t>{9) = 1 

^ ^ \2 2-Ki) \ 2 N 2-Ki) 

The remaining integral in H8U|) may be performed (similar as in |17]) with the 
result^ 

Similar as in |17] one can prove the eigenvalue equation 

dxT^\x)-Y,pf\ \9i,...,9n)t = ^ 
i=i J 

for arbitrary states. 

The fields ^^(x): Because the Bethe ansatz yields highest weight states we 
obtain the matrix elements of the spinor field V'l^;) = ^ipii^)- The p-function for 
the local operator il:^^\x) is (see also (HI]) 

P^'^\9_,z^=.^V±\(±^.-[l-^)±9^ ■ 

For example the 1-particle form factor is 

(0|V'(^)(0 ) \9)o. = Sal e^Hi-^)^ . 

^In |52II35| this result has been obtained using Jackson type integrals. 
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The last two formulae are consistent with the proposal of Swieca et al. |45l I38| 
that the statistics of the fundamental particles in the chiral SU (N) Gross- 
Neveu model should be o" = exp(27rzs), where s = ^ (l — is the spin (see 
also ()71() ) . For the n = N + 1 particle form factor there are again ni = N — I 
integrations in the ^-th level of the off-shell Bethe ansatz and the SU{N) weights 
are = (2, 1, . . . , 1, 1). Similar as above one obtains the 2-particle and 1-bound 
state form factor 

Km = N^e^-^i'-^)^'^ ^ - z)m - z)m - 

We were not able to perform this integration. In jSHl we will discuss the 1/A^ 
expansion. There we will also discuss the physical interpretation of the results 
for the chiral Gross-Neveu model. 
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Appendix 

A A Lemma 

Lemma 6 Let v^---" g V^ "^ be a highest weight vector, i.e. = 0. 

Then v^---'^ vanishes if the components tj^i^a-.-on Jq^. q,-^ = x vanish. 

Proof. By definition E acts on the basis vectors e^'"" = e^^ Cq-j (g • • • ea„ 

as 

n 

-Ei...„e^--" = ^ (g) • • • (g) Eca, (g • • • (g eQ,„ 
1=1 

Eca = e.a-1 where cq = . 
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Therefore we may write 

a,ai >1 
n 

0,Q!1>1 i=l 

o,«i=2 «,oi>2 

for some vector w^'"^. Because the second term w^'"^ is orthogonal to the first 
one ah components t>^°2 - on vanish. Iterating this procedure proves the claim. 



B Proof of Theorem [5] 

The proof of the main theorem of this article for SU{N) is a straightforward 
extension of the one for SU{3) above. 

Proof. The form factor equations (i) and (ii) may be proved quite analo- 
gously as for SU{3). The proof of (iii) is similar to that for the Z{N) model in 
[3l| . We use the short hand notations 

^=(01,...,e„), ^= (01,...,^^), ^= (0jv+l,...,^n), 

a = (ai, . . . , a„), a = (qi, . . . , an), a = (oat+i, . . . , a„), 

Z = {zi,..., Zm), Z= {Zi,..., ZN-l), Z = {zn, • • • , Zm), 

That F^' ^{6) given by (UHl), (El, (El and jZl satisfies (iii) in the form of ^ 
is equivalent to that i^i...n(^) satisfies 

n N 

Res . . . Res = Co TT TT ei...ArKAr+l...n(0 

eiv_lJV=«»? di2=lV • ^r , 1 o 

i=N+l j=2 

X {1 — aNSNn ■ ■ ■ Snn+i) (81) 

co = 2iV2''''TllF-^''~^\jirj) 

where the relation of F{z) and (/)(z) given by ()73() has been used. The residues 
of Ki ri{9) consists again of three terms 

Res . ... Res ^i...n(i) = Ml + 41 + 41 

because — 1 of the z integration contours will be "pinched" at three points. 
Again due to symmetry it is sufficient to determine the contribution from the 
zi, . . . , ZAT-i-integrations and multiply the result by m . . . (m — + 2). The 
pinching points are 
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(1) zi = 02, . . . , ZN-i = On, 

(2) zi = 01, . . . , ZN-i = On-i, 

(3) zi = 02 - ir], . . . , ZN-i = On - ir], 

The contribution of (1) is given by — 1 integrations along small circles 
around zi = 62, Z2 = O3, . . . , zn-i = On (see figure P). The S-matrices 5(^2 — 
zi), . . . , S{On — zn-i) yield the permutation operator 5(0) = P. Therefore for 
O12, ■ ■ ■ , 0n-2N-i,0n-in iv 



(n {0, z„) • • • Cf'^iO, zn)&''~' (0, On)--- {0, 

Pi Pm 




j=N 



Zj) ( Sin ' 



5i2 :' (n&^-{0,Zm)---C^''{0, 



ZN, 



- (82) 



It has been used that due to the SU{N) ice rule only the amplitude b{-) 
contributes to the S-matrices S{Oi — Zj) and a(-) to the S-matrices 5(^2 — 
Zj), . . . , S{ON — Zj), S{Oi — zi), . . . , S{0i—ZN-2) after having applied Yang-Baxter 
relations. One observes that the product of S-matrices in ()74p together with 
the one in ()82)1 yields the total — 1 S-matrix 



Sf-\z)iSiN...Sl2)-^ 



s 



TV.. .21 

a 



(0) 



for which the residue formula H69|) applies 



Res . 

di2=iv 



Res S, 



N...21 
ai...ajv 



{N - l)\{ir]) 



N-l^N...21. 



■^ai...ojv 



We combine (|82|) with the function Lp{z) with the property (|74|) for (/ = 1) 

and the scalar functions h and p and after having performed the remaining 
Zj-integrations we obtain 



r: 



(1) 

l...n 



n N 

« n n 

i=N+l j=2 



j) £l...N{0)KN+l...n{0) . 



We have used the following equations b{Oi — z)<f>{Oi — z) = —(j){z — O2) (for 
O12 = iv) together with the definition of t{z), the relation (iii') of H78() 
for the p-function p{0,02, . . . ,On,z) = (-1)™"^+^ p(il) (for Ou = --- = 
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9n~ 



IN 



irj), the relation (|79() for the normahzation constants NnCi{N 



l)l{iri)^^^ Y\f=i^ 4'U'i-v) = ^n-NCo and the recursion relation ((75|) for the con- 
stants q with the solution (|8H|) . 

The remaining contribution to H81|) is due to i?2 and R3 



r(2) , p(3) 



n Af 

-0 n n 

i=N+lj=2 



) £l...NKN+l...n{9)crNSNn ■ ■ ■ SnN 



+1 



It is convenient to shift the particle with momentum 9 i\f to the right by applying 
S-matrices and write 



n N 



+ 5^+1^ ...SnN + CO n n '^(%) ^l-NKN+l...n{9)aN = 



i=N+lj=2 



(83) 

where the components of this co- vector are now denoted by Wai...ojv-ijv+i---an«jv 
Note that because of (i) (fTT|) 



Kl,,,n{9)SN+lN ■ ■ ■ SnN = «(6'iAf )-^l...Ar+l...n7V (^1 , • • • , ^A^+l, ■ ■ ■ ,9n, 9^) ■ 

i=N+l 

Again due to Lemma El in appendix El it is sufficient to prove equation (|83j) 
only for = 1 since the left hand side of this equation is a highest weight 
co-vector. This is because Bethe ansatz vectors, in particular also off-shell 
Bethe ansatz vectors are of highest weight (see Therefore we consider 

(2) 

this equation for the components with = 1 only. The contribution of R\ ^ 
is given by the zi, . . . , ZAT-i-integrations along the small circles around zi = 
9i,. . . , ZN-i = 9n~i (see again figure^]). Now 5(^i — zi)i, . . . , S{9n~-i — zn~i) 
yield permutation operators P and the co-vector part of this contribution for 
a^v = 1 is 



(q C^™ {9, z„) • • • C^^-^ {9, 9n-i) {9, ^i)Pjv (1)) 

Pi Pn-I Pn 



•ajV-lJV+i---Onajv 



11 11 



^N-1 



i 








I3n- 


■.0 







-1 
-1 
-1 

9m 



where -P/v(l) projects onto the components with = 1- We have used the 
fact that because of the SU{N) ice rule the amplitude a(-) only contributes to 
the S-matrices S{9i — zj), S{92 — zj), S{9iy — zj), S{9i — zi) after having applied 
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Yang-Baxter relations. We use (l){6) = —b{9 + 2TTi)(j){9 + 2-Ki) to replace for 
i = l,...,iV- 1 and 1 

therefore using again (|(i9|) we obtain 

Res • • • Res 4>{0m) ■ ■ ■ 4>{0nn-i)S^:% ^(^i, • • • ,On-i) 

X Res • • • Res S^^^^i^ ^ (On + 27ri, , • • • , ^^v-i) 

6'jv-ijv="7 6i2=iV 

= -4>{iv) ■ ■ ■ HiN - l)iv){N - mivf-'ea,-a^,,i ■ 

We combine (|84|1 with the function Lp{z) with the property (|74|) (for / = 1) 

and the scalar functions h and p and after having performed the remaining 
Zj-integrations we obtain 

n N 

nfl^SN+lN ■ ■ ■ SnNPNil) = -CO J] 11 '^(%) £l...7v(l)i^7V+l...n(l)cTjV^Jv(l) • 

i=N+l j=2 

We have used the following equations: a{6iiy)4>{0ii) = 4>{0iN) (because of 
((T31) '). the definition of r(z), the relation (iii') of ((75|) for the p-function 
p{9_^ 9i, . . . , 9n-i,z) = (jp{6,z) (for 612 = • • • = ^at-iat = iry), the relation (|79|) 
for the normalization constants NnCi{N — l)l{iri)^^^ T[f=i <^(i*^) = ^n-NCo 
and the recursion relation (f75|) for the constants q. 

The contribution of pinching (3) is given by the zi,. . . , ZAr_i-integrations 
along the small circles around zi = 62 — ir], . . . , zj^-2 = Gn-i — zn-i = 
On — (see again figure^). As for = 3 the S-matrix S^p{dN — zn-i) yields 

r(^°") and the residue of the Bethe ansatz state vanishes, 

(pa) afi ' 



Res (nC{9,z^)---C{9,zi)PN{l) 

because is antisymmetric with respect to a, 13. Therefore equation ()83() is 
proved for ajy = 1 and because of Lemma El also in general. ■ 

C The higher level Bethe ansatz 

Proof of lemma |3J For the higher level functions L^'"* (z) one may verify the 

equations (i)*-'^ and (11)^'-' quite analogously to the corresponding ones for the 
main theorem (e.g. for = 3). 

We prove (iii)^') by induction and assume: 

k N-l-1 

if^)(2x)^Q+i5f-'+2(^) Jl Jl 4>{u^,)lI~''\u) (85) 

i=N-l j=2 
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for ui2, . . . , UAr-i-2,Ar-2 ^ "^V- t^e integral representation ((7^ of L^j^ {z) 
there are pinchings at ui = Z2, • • • , un-1-i = zn-i if 212, ... , zn-i-i,n-i ifl- 
Therefore in {z_,u) the S-matrices S{z2 — ui), . . . , S{zj\[-i — uat-z-i) yield 
the permutation operator P and we have to consider 

^f-{z,u) = (ow • • • c-w^^-'-d, z^.i) ■ ■ ■ c-w^Hi, Z2))^ 

1+1 1+1 1+1 



Ik 

71 























Uk 






















UN-1-1 








ZN-l 


Ul = Z2 


Z-m 



1+1 

1+1 
1+1 



(3n 



{Sl2{zi2) . . .SiN^l{ziN~l))' "^"^ ' ' TT b{zi-Uj)^^^)-{z,u) 



where ^ = (/^AT-z+i, ... ,/3m), 1 = (zn-i+i, ■ ■ ■ , Zm), u = {un-i, ■ ■ ■ ,Uk). We 
may write for « = (m, . . . , un-i-i) = (z2, • • • , zn-i) 



s 



N...I+1, 

'/3 



with the notation z_= (zi, . . . , zj\[^i). Therefore using the assumption we 
obtain when Z12, ■ ■ ■ , z^^i^i^^^i irj 



dui 



kl Jz2 ^ . ^ 



dUN-l-l f duN-i 



m N~l 



h{z_, u) 



X L('+^)(n)l> 



^^Hz,u)=c,si-^^Hz) n n^(-^.-)4'^(i) 

- j=Ar_;+i j=2 - 



where k = k — N + I + 1. The following formulae have been used: b(zi 



{zi-Uj)l{i=2 -rizi-Uj) 



= — 1, the relation (iii')^^-* of 
JZZI) for the p-function p'^^ki^, Z2, . . . zn-i,u) = (-1)'=~^"'+Vjj^\(l, n) (if zu = 
_i N-i ~ recursion relation q = q+i H^ITl'^^ The 



solution to this recursion relation with cn-i = 1 is 



(86) 
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